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Abstract 

In this series of lectures we present the functional integral method 
for studying the superconducting pairing of quarks with the forma- 
tion of the diquarks as well as the quark-antiquark pairing in dense 
QCD. The dynamical equations for the superconducting order pa- 
rameters are the nonlinear integral equations for the composite quan- 
tum fields describing the quark-quark or quark-antiquark systems. 
These composite fields are the bi-local fields if the pairing is gener- 
ated by the gluon exchange while for the instanton induced pairing 
interactions they are the local ones. The expressions of the free 
energy densities are derived. The binding of three quarks is also 
discussed. 

1 Introduction 

The superconducting pairing of quarks due to the gluon exchange in QCD with the 
formation of the diquark condensate was proposed by Barrens' 1 ' and Frautschi' 2 ' 
since more than two decades and then studied by Bailin and Love' 3 ' , Donoglue and 
Sateesh' 4 ', Iwasaki and Iwado' 5 '. Recently, in a series of papers by Alford, Rajagopal 
and Wilczek' 6 !, Schafer and Wilczek' 7 !, Rapp, Schaier, Shuryak and Velkovsky' 8 ! , 
Evans, Hsu and Schwetz^, Son' 10 !, Carter and Diakonov' 11 ! and others' 12-21 ' there 
arose a new interest to the existence of the diquark Bose condensate in the QCD 
dense matter-the color superconductivity. The connection between the color super- 
conductivity and the chiral phase transition in QCD was studied by Berges and 
Rajagopal' 22 ', Harada and Shibata' 23 '. There exists also the spontaneous parity 
violation, as it was shown by Pisarski and Rischke' 24 ' . 

For the study of many-body systems of relativistic particles with the internal 
degrees of freedom as well as with the virtual creation and annihilation of the 
particle-antiparticle pairs the functional integral technique is a powerful mathemat- 
ical tool. This method was applied to the study of the color superconductivity as 
well as the quark-antiquark pairing' 25-28 ! , 

In this series of lectures we present the basics of the functional integral method 
for the study of the color superconductivity in QCD at finite temperature and 
density (the so-called dense QCD or thermal QCD). As the physical origin of the 
superconducting pairing of quarks we consider two different commonly discussed 
mechanisms: the direct local four-fermion interactions of the quark field and the 
quark-quark non-local interaction due to the gluon exchange. The direct four-quark 
interaction might be induced by the instantons' 11,12 '. 

We follow the general method of the functional integral approach in the theory 
of superconductivity I 29 ' 30 ] . We start from the expressions of the partition functions 
of the systems of interacting quarks and antiquarks with some local or non-local 
quark-quark interactions. Then we introduce the local or bi-local composite fields 
describing the diquarks or quark-antiquark pairs, establish the effective action of 
these composite fields and derive their field equations. The order parameters of 
the ground state of the diquark or quark-antiquark condensate are the expectation 

1 Lectures at the Vlth International School in Theoretical Physics, Vung Tau 27 December 
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values of the composite fields in the corresponding state of the system. Due to the 
translational invariance of the ground state these expectation values are coordinate- 
independent (for the local fields) or depend only on the difference of the coordinates 
(for the bi- local ones) . The equations for the order parameters are the special cases 
of the field equations. The expressions of the free energy density in the corresponding 
phases are also derived. 

We work in the imaginary time formalism and write briefly 





x = (x, r) , J dx = j dr J dx, 



A 

k is the Boltzmann constant and T is the temperature. We denote tp A , ip the 
quark field and its Dirac conjugate, where A = (a,a,i) is the set consisting of the 
Dirac spinor index a = 1, 2, 3, 4, the flavor index i = 1,2, 3, ....Nf and the color one 
a = 1,2,3...N C . The internal symmetry groups are assumed to be SU (Nf)j, and 
SU (N c ) c . The partition function of the system of free quarks and antiquarks with 
the chemical potential [i and at the temperature T can be expressed in the form of 
the functional integral 



Z = J [Dip] [Dip\ exp 



where 



dxtp A (x)L B A i) B (x) j 



/i + 7V + M 



(1) 



(2) 



and M is the bare quark mass. 

Introduce the generating functional 



Z [ V ,rj] = J [£>V] [D$\ exp j- J dxyj A (x) L A tp B (x) 



exp 



dx 



V (x) i> A (x) + ip {x) r) A (x) 



(3) 



with anti-commuting parameters i] A (x) and r] A (x). The 2n-point Green functions 
are expressed in terms of the functional derivatives of Z [77,77] at the special value 
Va (x) =r} A (x)=0: 



g a\','.'.'.'.aZ i x i,x2, -,x n ;yi,y2, -,y n ) = 



= (-1) 



{xi)ip A2 {X2) —i>A n (x n )ip Bn (y n ) ...ip B2 {y2)4> Bl (yi)}) 

5 2n Z [77,77] 



Z Sr) 1 (x!)...5r) " (x n )5r) Bn (y n ) ...Sr} Bl (j/i) 



\r}— r\— • 



(4) 



In particular 



G B A (x-y) = G B A {x-y) = (r{i> A {x)^ B {y))) 
1 S 2 Z [t,,t}] 

Z SrJ A (x)Sr lB (y) 1 ^ - (b) 

Denote 

S% (x-y) = 5 b a 5lSi (x-y) = 5 b j{S^ (x - y,r - a) , (6) 
Si (x - y,r - a) = S afj (x - y,r - a) 



L B A S c B {x-y) = 5 c A 5(x-y), (7) 
z = (x,t), y=(y,cr), <5(x-y) = 5 (x - y) 5 (t - a) , 



7 4 (J^ - f^j + 7V + M S (il (x - y,r - cr) = <5 Q7 <5 (x - y) 5 (r - cr) . 

(8) 

Shifting the functional integration variables 

♦.(.) - *.W + /*sg(.-»), D W, 

(x) -» ^ A (a:) + /" dzrf (z) S£ (z - x) 

in the r.h.s. of the formula (1), we derive the explicit expression of the generating 
functional (3) 



Zq [v, v] = Z exp jy dx J dyr\ A (x) S% (x - y) rj B . 



(9) 



Substituting this expression into the r.h.s. of the formula (5), we obtain the two- 
point Green function 

G B A (x-y) = (T [i, A (x) ^ (y)}) = S% (x - y) . (10) 

Similarly, from the formulae (4) and (9) it follows the Wick theorem for the 2n-point 
Green functions of the free fermionic fields. 

Introduce the Fourier transform S a p (p) = S a p (p, e m ) of S a p (x) = S a p (x, r) : 



Sap (X, T) = 1 £ 6 " mt ^3 y d P ^ S «P (P. £m) 



(11) 



(2m+l)|, 



m being the integers m = 0, ±1, ±2... From the equation (8) it follows that 

^ (p)= 1¥TmJ q/3 = p2 + M < 12 > 

with the notations 

P= (e m + «m)74 + 7P, P 2 = (£m +«/") 2 + P 2 - (13) 

In the calculations we shall use also the expression of S a p (—p) = S a p (— p, — e m ) . 
For the convenience we write it in the form 

5 < ^ ( 1 ^ ^ + M U ,,,, 

Safi { - p) = \^WTVi) aP = p2 + M ■ (14) 

where 

= (e m - in) 7 4 + 7p, p' 2 = (e m - z^) 2 + p 2 . (15) 

To study the color superconductivity we consider the quark-quark pairing in 
QCD with the formation of the diquark condensate. Then we investigate the quark- 
antiquark pairing. The corresponding composite fields are the meson ones. We 
discuss also the possibility to extend our reasonings to the study of the binding of 
three auarks. The comnosite oarticles in this case are the barvons. 



2 Quark- Quark Pairing 



For the simplicity in writing the formulae we begin our study by considering the 
superconducting quark-quark pairing due to some direct four-fermion coupling of 
quarks with the interaction Lagrangian 



Lint = ^ A {x)^ C \x)Vg£ip D {x)ip B (x), 

\rBD _ \rDB _ \rBD _ \rDB 
V AC — V CA — ~ V CA — ~ V AC ■ 



(16) 



The partition function of the system equals 



exp < — / dxip (x)L A ip B (x) 



-rA. x — c 



dx^\x)^(x)V^iP D (x)iP B (x) } . 



-=AC 



Introduce the antisymmetric bi-spinor local fields &ca (x) , $ (x) 

$ (x) = 

and the functional integral 



$ CA (x) = -$ AC (x) , $ac (x) = -<$>ca (x) , 



(17) 



(18) 



Zq = J [D$] [DW] exp J dx<S> A ° 



{x)Vg2{x-y)* DB {x)\. (19) 
By shifting the functional integration variables 

$db(x) -» $db (x) + tp D (x)ip B (x) , 
(x) — > ¥ AC (or) + il> A (x) ip C (x) , 
we establish the Hubbard-Stratonovich transformation 



CXP {^ / ^ ^ ^ V ca^d (z) ^b ( x ) 

± J [£>*] [£>$]exp|-i J dx<p AC {x)V£2* DB (x)} 



(20) 



where 



-BD 



-AC 



A— (x) = *™ (x) Vgf, A CA (x) = VgZ* DB (x) , 

A° B (x) = -A BD (x) , A AC (x) = -A CA (x) , 
and rewrite the partition function (17) in the form 

Z= ^fJ m] ^ 6XP ^° ff [$ '^ } 
with the effective action 

S cS [<&,$]= -\ J dx$ BD (x) Vgg*cA (*) + W [A, A] , 
exp {W [A, A]} = 

= i + £r( 2 ") [a, a] 

n=l 

= /t [exp (-1 / dx fe^ (x) ^ (x) A CA (x) + A B ° (x)^ D (x) V B ( 



(21) 
(22) 

(23) 
(24) 



(25) 



Calculations give 

[A, A] = W {2) [A, A] = (26) 
= \j dx ^j dx 2 A AlCl {xi) S% {xi - x 2 ) A C . 2 A 2 (x 2 ) Slf 2 (x 2 - xi) , 

[A, A] = i (^) [A,A]) 2 + ^ 4 ) [A, A] , 



(27) 



W (4) [A, A] = ~\J d Xl ... J dx 4 A AlCl (xi) S% (xi - x 2 ) A C2 a 2 (a*) S™ 2 (a* - a*) 



A A3 ° 3 (x 3 ) Sg* (x 3 - Xi ) A C4 a, (a*) SXt 4 (a* - a=i) 
r< 6 > [A, A] = 1 (V 2 > [A, A])' + [A, A] If* 4 ' [A, A] + W* 6 ) [A, A] 



(28) 



(29) 



W (6) [A, A] - 1 /" dan... /" cfeeA- 4101 (an) Sg 2 (a* - x 2 ) A C2 .4 2 (ar 2 ) S™ 2 (a* - x 3 ) 
... A- 4505 (x 5 ) (x 5 - x 6 ) A C6 A e (x 6 ) Slf e (x e - xi) , (30) 



with 

We have then 



Sl A (jci - ar 2 ) = S| (x 2 - xi) . 



W [A, A] = w(2n) [ A ; • ( 31 ) 

n=l 

From the variational principle for the effective action we derive the field equation 



1 A M yBn^ SW^ [A, A] 

-A C A (X) = V CA 2^ BD 

1 „ = i 6A (a;) 



(32) 



Using the explicit expressions of W^ 2 ™) [A, A] , we obtain 

A CA (x) = V^A I y dx 2 S c D 2 {x-x 2 ) A C2A2 (x 2 ) S T B M (x 2 - x) 

- j da*..- y d^S'g 2 - x 2 ) A C2 a 2 (x 2 ) Sjf 2 {x 2 - x 3 ) 

A A3Cs (x 3 ) (x 3 - x A ) A Ci Ai {xa) Sg Ai (xa - x) 
+ y da*... y <fe 6> s£ 2 (a: - a*) A C2j42 (ar 2 ) Sjf 2 (ar 2 - a*) (33) 



A 3 ° 3 (x 3 ) Sg< (ar 3 - xa) ... A Ce .4 6 (a*) (x 6 - x) 

} 



In the special class of the constant solutions 

A CA (x) = A CA = Vgf* DB , A B ° (x) = A BD = ^ AC V^, (34) 
we have the extended BCS gap equation 



(27T)' 

S(p,e m )AS T (-p,-e m ) 



l + A5( Pl e ra )AST(-p- £m ) 



,(35) 



where S(p,e m ), S T (— p, — e m ) , A and A are the matrices with the elements 

$A (Pi £ m) ! 



AC 



-BD 



S\ B (-p,-e m ), A c a, and A 

At the values of the superconducting order parameters Ac a and A^^ satisfying 
the extended BCS gap equation the effective action equals 

5eff[$,*] = (\-\)J dx i- J dx 4 Tr[AS(x 1 -x 2 )AS T (x 2 -x 3 ) 

A S(x 3 ~x i )AS T (x i -x 1 )] (36) 

~ G~ o") / dai1 "" / dx 6Tr[A5(a; 1 - a ; 2 )A5 T (x 2 - x 3 ) 
...AS(x 5 - x e )AS T (x 6 - Xl )] 



Denote F [x; A] the free energy density of the condensate. The effective action is 
expressed in terms of this free energy density in the following manner 



S off [$,$] = -/3 y"dxF[x;A]. 
Comparing (36) with (37), we obtain 



(37) 



F[x; A] = F [A] = / 



dp 



Tr | Q - ^ [A 5 (p, e m ) A S T (-p, -e m ) 



1 f 

2 ~ 6 



•f f 
- [ 2 - 8 



AS(p,£ m )AS T (-p,-£ m ) 



AS(p,gA5 T (-p,-£ m ) 



(38) 



Summing up the infinite series, we write the r.h.s of the formula (38) in the compact 
form 



F[x;A] = ^V — 
$ ™ (2*) 



S (p, e m ) AS T ( p, —e m ) A 



1 



1 + S (p, e m ) AS T (-p, -£ m ) A 
i 

c?a- 



(39) 



f 



o 



1 + aS (p, e m ) AS T (-p, -e m ) A 



Let us discuss the general form of the superconducting order parameters. Con- 
sider first the constants Aac- We have 



A A c = A (aia)(ckl) = (75C) AL )(cfe) + (C) A 



/Q7 (oi)(cfc) ' 



(40) 



where A?^,^ are the scalar constants while the Af^,^ are the pseudoscalar 
ones. If the parity is conserved, then all pseudoscalar constants Af^w^N must be 
zero. The existence of non- vanishing pseudoscalar constants A^\^ c ^ would signify 
the spontaneous breaking of the parity conservation. Because of the condition (22) 
the constants A^^ and A^^ cfe ^ must have the property 



A 5 



A 5 



A P 



^(ai)(cfc)- 



(41) 



This means, in particular, that if they are symmetric (antisymmetric) with respect 
to the flavor indices i and j, they must be also symmetric (antisymmetric) with 



For the study of the quark-quark pairing due to the gluon exchange we start 
from the partition function in the form 



where 



J [Lhl>][D1>] ex P j-/ dx^ A {x)L B A i, B {x) 
.expji Jdxf dy^ A (x)^ C (y)V^ (x - y)^ D (y)^ B (x)Y (42) 



D 1 



Vca (x y) 27r2 ^UC-'M wc-'/c (x-y) 



Ma - Jf d a d c 



,(43) 



A/ are the Gell-Mann matrices of the color symmetry group. In order to describe the 
diquark systems we introduce some composite bi- local bi-spinor fields <&bd {x,y) , 



—AC 



$ (x,y) obeying the Fermi-Dirac statistics 

<f> DB (y, x) = -<S>bd (x, y) , (y, x) = -$ AC (x, y) , 



(44) 



and the functional integral over these bosonic fields 

Zt = J [£>$] [£>*] cxp |-1 y dx J dy$ AC (x, y) Vgf (x - y) <S> DB (y, x) J . 

(45) 

By means the shift of the functional integration variables 



®db (y,x) 



$db (y,x) + ip D {y)ip B (x) 



$ AC (x,y) -> $ AC (x, y) + tp A (x) ^ (y) , 
we can establish the Hubbard-Stratonovich transformation 



(46) 



^ y [£><&] [£>$]exp|-l I dx y dy$ AC fry)^ (x-y)$ DB (»,*)} 



cxp \ -i / r/.r / tin 



A (x, y) ^ D (y) ip B (x) + ip {x)ip (y)A C A(y,x) 



}■ 



A C A (y, z) = (x - y) $ DB (y, x) 



-BD 



^AC 



(47) 



A AC (x,y) = -A C A(»,a;), A DS (y, x) = -A B ° (x, y) , (48) 



A {x,y) = * (x,y)V*°{x-y), 
:), A 
and rewrite the partition function 

z = ~zj I m] cxp {A I dx I dy ^ AC {x ' y) {x ~ y) ® db (y ' x) ] 



j [Dip] [Dip] exp|-y 



dxip (x) L^iPb (x) 



(49) 



in the form (23) with the effective action 
S cS [<&,$]= -\ J dx J dy$ AC (x, y) Vq® (x - y) $ flB (y, x) + W [A, A] , 



} 



= T 



cxp{W [A, A]} = 

OO 

1 + ^ r (2n) 

1 



n=l 



exp 



dx / c?y 



A {x,y)4> D (y)il> B (x)+i> (x)i> {y)A C A(y,x) 



Calculations give 

r^) [a, a] = w* 2 > [a, a] = i y da* y di/i y dx 2 y ^ (52) 

A AlCl (xi,2/i) Sg 2 - j/ 2 ) A CaA2 (ys,^) S T A ^ (x 2 - Xl ) , 
S T A ^(x 2 -x 1 ) = Sil(x 1 -x 2 ), (53) 
r W [A, A] = 1 (V< 2 > [A, A] ) 2 + ^< 4 ) [A, A] , (54) 

1 



[A, A] = 



ctei / dyi / dx 2 / dy 2 / dx 3 / rfy 3 / dx 4 / dj/4 



(xi, yi ) S^ 2 (2/i - 2/2) &c 2 a 2 (2/2, £2) S™ 2 (x 2 - x 3 ) (55) 



(2:3,2/3) Sg 4 (^3 - 2/4) A C4 a 4 (2/4,2:4) S 1 ^ 4 (x 4 - Xi) , 



TA 4 



[A, A] = i (V 2 ) [A, A])' + W< 2 ) [A, A] WW [A, A] + W< 6 > [A, A] , 



(56) 



[a, a] = i y dxi y y dx 2 y ^ ... y y ^ y ^ 6 y 



^2/6 



-A1C1 



(xi, 2/1) S^l (2/1 - 2/2) A C2 a 2 (j/ 2 , x 2 ) S^ 2 (x 2 - x 3 ) 



-T-A 3 C 3 



(%3, 2/3) Sg 4 (2/3 - 2/4) A C4 a 4 (2/4, £4) 51 j^ 4 (^4 - 2:5) (57) 



T A 4 



-A5C5 



(2:5, 2/5) S%1 (2/5 - 2/6 ) A C6 a 6 (ye, x 6 ) S^ 6 (x 6 - xi) , 



Ti 6 



r( g )[A,A] = i(ff( 2 )[A,A]) 4 + i(^[A,A])V( 4 »[A,A] 



(58) 



~ (W< 4 » [A, A])' + VK< 2 ) [A, A] W^ 6 > [A, A] + ^ [A, A] , 



H^ (8) [A, A] = y^X! y dy 1 J dx 2 J dy 2 ... j dx 7 j dy 7 J dx s J dy 8 



-A1C1 



A" 1 " 1 (x l7 2/1) S% (2/1 - 2/2) A C2 A 2 (2/2, 2:2) S^ 2 0»* - 2:3) 



TA 2 



A" 3 "" 3 (x 3 , 2/3) S# (2/3 - 2/4) A C4A4 (2/4, x 4 ) S™* (x 4 - x 5 ) (59) 



-*-A 5 C 5 



(2:5, 2/5) S^l (2/5 - 2/6 ) A C( , A6 (j/6, x 6 ) S^ 6 (x 6 - x 7 ) 



A A7C7 (x 7 , 2/7) Sc* (2/7 - 2/8 ) A C8 a 8 (2/g, x 8 ) S^ 8 (x 8 - xi) , 



It is easy to verify that 



W [A, A] = J2 w(2n) [ A > 2] ■ 



(60) 



From the variational principle for the effective action we derive the field equation 

l -A CA (y,x) = Vg? (x-y)J2 - BD [ ' J - (61) 
z n= i 5 A (x,y) 

It has the explicit form 

A CA (y, x) = VgZ {x~y){j dx 2 J dy 2 S% 2 (y - y 2 ) Ac 2 a 2 (y2, x 2 ) S t b m (x 2 - x) 

- J dx 2 J dy 2 ... J dx 4 J dy A S^ 2 (y - y 2 ) A C2 a 2 (y 2 ,x 2 ) S^f 2 (x 2 - x 3 ) 
A A3 ° 3 (x 3 , y 3 ) (y 3 - y 4 ) A CiA4 (y 4 , a*) S T B A * (x 4 - x) 



+ J dX2 J dm '" J dX6 J dy6S ° 2 ^y~V2)A C2 A 2 (y 2 ,x 2 )S Aa 2 (x 2 -x 3 ) 
A AsC3 (x 3 , y 3 ) (y 3 - y 4 ) A CiAi (y 4 , x 4 ) S^ Ai (x 4 - x 5 ) 
A A5CS (a*, y 5 ) Sgj (y 5 - y 6 ) A CfiA<i (y 6 , x 6 ) S^ 6 (a* - x) 

- j ' dx 2 J dy 2 ... J dx 8 J dy 8 S% 2 (y - y 2 ) A C2 a 2 (y 2 ,x 2 )S^ 2 (x 2 - x 3 ) 

... A ArC7 (x 7 , y 7 ) (y 7 - y 8 ) A C8 a 8 {Vs, x 8 ) S£ a * (x 8 - x) 
+ ... }. (62) 

Considering the solutions of this equation in the special class of functions depending 
only on the difference of the coordinates 

Aca (y, x) = A c a (y - x) , A BD (x, y) = A BD (x - y) , (63) 
performing the Fourier transformations 

A c a (y - x, a - t) = - V e^iv-T) _L_ f dpe^-^A C A (p, e m ) , 

P m ( 27r ) J 

A BD (x-y,T-a) = W^^Tr^ /dpe ip(x - y) 5 SD (p, £ro ), (64) 
V^f(x-y,r-<7) = ^ e " ra( '" T) 7^/ d P e<P(y " X) ^ B (P,Wn,), 

TT IT 

£=(2m + l)-, uj m = 2m-, 
and introducing matrices A (p, e m ), A(p,e m ), S(p,e m ), S T (— p, — e m ) with the 

„ ^_ BD „ „ 

elements A C a (p, s m ), A (p,e m ), S A (p,s m ), S], A (-p, -e m ), we derive the 
extended BCS gap equation for the superconducting quark-quark pairing in QCD 

A CA {p,s m ) = \Y\— \s [ dqVcA {p-q,e m -e n ) 

\S (q, e n ) A (q, e n ) S T (-q, -e„) 

.(65) 



f 



1 + A (q, e n ) S (q, £„) A (q, e n ) S T (-q, -£„) 



At the values of the fields Aca (y, a;) and A B ° (x, y) satisfying the equation (62) 
the effective action (50) equals 

S eff [$, $] = W[A,A\- [dxf dyA AC (x, y) ^Jf ' ^ 



= J dxi J d Vi--- J dx ± J dy 4 Tr[A(x- L ,y 1 )S(y 1 -y 2 )A(y 2 ,x 2 ) 



+ 



S T (x 2 - x 3 ) A (x 3 ,y 3 ) S (y 3 - y 4 ) A (y 4 , x 4 ) S T (x 4 - zi)] 

J dx i j dyi- J dx 6 J dy 6 Tr [A (xi,yi) S (j/i - y 2 ) A (y 2 , x 2 ) 

S T (x 2 - x 3 ) ... A (x 5 ,y 5 ) S (y 5 - ye) A (y e ,x e ) S T (x 6 - xi)] 

J dxi J ' d Vi- J dx » J dy 8 Tr[A(x 1 ,y 1 )S(y 1 -y 2 )A(y 2 ,x 2 ) 

S T (x 2 - x 3 ) ... A (x 7 , y 7 ) S {y 7 - y$) A (y 8 , x 8 ) S T {x 8 - x L )\ 
+ ..., (66) 

where A (y, x), A (y, x) , S (x — y) , S T (x — y) denote the matrices with the elements 



-AC 



Aca (y, x), A {x, y), S*l (x — y), Sg (x — y). In the case of the special class (63) 

of the solutions Aca (y — x), A {x — y) the effective action is expressed in terms 
of the free energy by means of the formula (37). Comparing the expressions (37) 
and (66) and performing the Fourier transformation, we obtain 



Tr iU~i) F (p ' £m) (p ' £m) ^ (p ' £m) ^ T ( p ' ~ Em) 



1 1 

2 ~ 6 



'I 1 
'28 



A (p, £ m ) S (p, S m ) A (p, E m ) S T (-p, -£ m ) 



A (p, E rn ) S (p, £ m ) A (p, E m ) S T (-p, -£ m ) (67) 



This series can be also written in the compact form of the integral representation 

S (p, E m ) A (p, E m ) S T ( p, -E m ) A (p, £ m ) 
1 



F [A] = \ Y ( dpi Tr 



1 + S (p, e m ) A (p, e m ) S T (-p, -£ m ) A (p, e m ) 
i 

— / da- 



1 



o 



1 + aS (p, £ m ) A (p, £ m ) 5" T (-p, -£ m ) A (p, £ m ) 



(68) 



In order to simplify the expression for Aac (p ; £m) we introduce 4- vector p u with 
the fourth component p 4 = Em and denote Aac (p, £m)by Aac {p)- The condition 
(48) is now written in the form 

a ba {-p) = -A AB (p) ■ (69) 
The bi-spinors Aab (p) are expressed in terms of the Dirac bi-spinors 

A AB (p) = A [aia)(bjp) (p) = (7 5 C) a( g Af oi)(M) (p) + (C) afj Af ai){bj) (p) 

(P) (70) 

+ {onvl b C) afj A^ v{ai){bj) (p) . 

From the conditions (69) it follows that under the interchange (ai) <-> (bj) and the 
inversion p —p the scalar, pseudoscalar and vector functions are invariant while 
the pseudovector and tensor functions change their sign 

A (bj)(ai) (~P) = A (ai)(bj) (P) ' A l(bj)(at) (~P) = A ^ (ai)(bj) b) ' 
A . , - -A A . . .. . (n\ A* . t- n \ — -A* . . 



In general, the existence of non-vanishing superconducting order parameters 
which are not the singlets of the color and/or flavor groups and lower the free energy 
would mean the spontaneous breaking of the color and/or flavor symmetries. For 
the systems with the isomorphic color and flavor groups SU (N) c and SU (N)^ 
there may exist the superconducting order parameters which are the irreducible 
spinor representations of the groups SU (N) c and SU (N)^ but the singlet of the 
"diagonal" SU (N) subgroup of the direct product SU (N) C ®SU (N) f . In this case 
we have the " color- flavor locking" . 

3 Quark- Antiquark Pairing. 

The direct four-fermion coupling of the quark fields with the interaction Lagrangian 
(16) or the non-local interaction of the quark fields with the effective action given 
in the r.h.s of the formula (42) are also the origins of the quark-antiquark pairing. 
In the case of the quark-antiquark pairing due to the direct four-fermion coupling 
of the quark field we use the interaction Lagrangian in the form 



-7 A 

2 



Lint = ^ (x) i> B {X) U%g 1p (x) 1}J D (x) , 



jjDB 
U CA 



jjBD _ ttDB _ jjDB 

~ u AC — U CA — u AC ' 



where instead of the constants V®a m (16) we use the new notations 

jjBD _ V BD 
° AC — V CA ■ 

The partition function of the system equals 

Z = J [Lh/>] [D$\ exp |- J dx^ A (x) L B A i, B (x) J 

dx^ (X) ^ B (X) UlB~f (X) Ipjj (x)\ . 



(72) 
(73) 



exp 

I z J 

Introducing the local hermitian fields $^ (x) and the functional integral 

Zt = J [£>$] exp |-1 J dx^i (x) U%g$% (x)| , 
we establish the Hubbard-Stratonovich transformation 

exp |^ J dxip A (x) ip B (x) U B E if (x) ^ D W | = 
= jwj l D *} cx v{-\J dx ^B (x) UlB* c D (*)} 
exp | - J dx~4) A (x) i) B (x) A| (x) | , 

and rewrite the partition function in the form 

with the effective action 

Seff m = ~\J dxH (x) U1S* C d (x) + W[A], 

exp {W{A]} = l + jr r (n) [A] = /t exp j- f dxA 

n=l \ L ^ ^ 



(74) 



(75) 



(76) 
(77) 

(78) 
(79) 



A (x) 1p (x) 1p B (x) 



Calculations give 

CO 

W [A] = w(n) i A ] > ( 81 ) 

n=l 

[A] - J dxA B A (x) S£ (0) , (82) 

[A] = -1 y | dx 2 A^ ( Xl ) S£ (xx x 2 ) A% (x 2 ) S% (x 2 - Xl ) , 

(83) 

[A] = l -Jd Xl ...J dx 3 A B A \ ( Xl ) s£\ ( Xl - x 2 ) A% (x 2 ) 



Sp 2 (x 2 - x 3 ) A B A l 0*0 S£l (x 3 - x,) , (84) 



W (n) [A] = ( ^ + [dxi... f dx n A B A \ (an) < 2 (si - a*) A*° (x 2 ) 

s£ («2 - a*) - Af » (s„) S£ (*„ - a*) . (85) 
From the variational principle 

SS eff (*) 



we derive the field equation 



= (86) 



AS(x) = Ug^^P- = Ug^(x,x), (S7) 
<5 A3 (a:) 



where (y, x) is the two-point Green function of the quark field in the presence of 
the pairing interaction 

(y, x) = S £(y-x)-j dx,S^ (y - Xl ) A*] ( Xl ) ( Xl - x) 

+ j " d Xl J dx 2 S£> (y- Xl ) A*l ( Xl ) s£\ ( Xl - x 2 ) Af 2 (x 2 ) S£ 2 (x 2 - x) 

- J dx x ... J dx 3 S^ (y - Xl ) A B A \ (xi) S^l (xx - x 2 ) A B A \ (x 2 ) (88) 
S£(x 2 -x 3 )A B l (x 3 )Sl(x 3 -x) + .... 
It satisfies the Schwinger-Dyson equation 

fa, x) = S£(y-x)- J dzS c B (y - z) Ag (z) s£ (z, x) . (89) 

In the special class of the constant solutions of the field equation (87) 

A^ (x) = Ag = const (90) 

(y, x) depends only on the coordinate difference 

Si(y,x)=8i(y-x). (91) 

For its Fourier transform we have then an algebraic equation. Denote A the matrix 
with the elements A^. From the equation (89) it follows that 

S (p, e m ) = S (p, e rn ) - S (p, e m ) AS (p, e m ) (92) 

and 

^ = ^ + A. (93) 



With the field satisfying the equation (87) the effective action equals 
S eff m = W[A]-\j dxA B A (x) ^± =( 1 -\)j dxA * & S * W 

+ {l~l)J dxl J dx2 J dx3 A ^ {xl) s £ {X1 - X2) A S {X2) 

S£l(x 2 -x s )Af 3 (x 3 )S£l(x 3 - Xl ) (94) 

dXl - J dXiA *> (Xl) S *' {Xl ~ X2) A % {X2) 
S£ (x 2 - x 3 ) ... A* (x 4 ) S£l (x 4 - Xl ) + ... . 

= ^K 1 -^) / dxA(x)S(0) 

+ {\ ~ ^) J dx 3 A (xi) S ( Xl - x 2 ) 

A (x 2 ) S (x 2 - x 3 ) A (x 3 ) S (x 3 - xi) 
~ (i ~ 2~) j dxi '- J dx 4 A(x 1 )S( Xl -x 2 )A(x 2 ) 

S(X 2 -X 3 )...A(X4)S(X4-X 1 ) + .... j, 

where A (x) is the matrix with elements A^ (x) . ft follows that in the case of 
the constant solutions (90) of the field equations (87) we have following formula 
determining the free energy density F [x; A] : 



F [x; A] — F [A] 



l--)AM(O) 



+ J d yJ dzA B A S c B {x-y)A°S*{y-z)A F E S*{z-x) 

JdyjdzJ dwA B A S C B (x - y) A B S% (y - z) 



= Tr 



1 1 

4 ~ 2 

A£ S$ (z - w) A%S£ (w-x) + 
I - I 1 AS(0) 



(95) 



+ (l-iyjdyj dzS(x-y)AS(y-z)AS(z-x) 

-(^-^jAjdyJdzJ dwS(x-y)AS(y-z) 
A S (z — w) A S (w — x) + ... 



(27T) 

1 l\ r ~ i 3 
3 2 J [AS(p,e m )_ - 

Summing up the infinite series, we obtain 



1 1 

4 ~ 2 



AS(p,£ ra ) 



F [ A ]^E(^/ d P Tr { A (P< £ ™) ^ 1 S«( P ,e m )^-is(p, £m ) j, 

(96) 

where S (p, e m ) satisfies the equation (93) and (p, e m ) is determined by a similar 
one with the replacement of A by £A: 



+ £A. 



(97) 



The order parameters A^ have the form 



(oi) 
P («) 



la (ai) 



(98) 



The non- vanishing order parameters A, \ ' lowering the free energy would mean 
the spontaneous parity conservation violation. If M = and the interaction La- 
grangian (104) is invariant under the chiral transformations, then the existence 
of non- vanishing order parameters A^J^ and/or A^j* lowering the free energy 



(ai) 



would signify the spontaneous breaking of the chiral invariance. If Af ( .f and/or 



arc not tnc sm gl e ts of the color and/or flavor group, then the color and/or 
flavor symmetries are spontaneously broken. For the system with the isomorphic 
color and flavor groups SU (N) c and SU (N)^ there may exist the superconduct- 
ing order parameters which are the irreducible spinor representations of the groups 
SU (N) c and SU (N) f but the singlet of the the "diagonal" SU (N) subgroup of the 
direct product SU (iV) c ® SU (N)^. In this case we have the "color-flavor locking". 

Now we consider the quark-antiquark pairing generated by the effective non- 
local interaction of the quark fields due to the gluon exchange. For this purpose we 
rewrite the partition function (42) in the appropriate form 



(ai) 



z = 



[D$\exp{-J 



dxip {x) L^ip B ( x ) 



\ [ ~ j ll >- j d y^ A (*) ^ B {y)u^S(x^y)^ C (y) i> D (x) } . (99) 



with the new notations 



U%?{x-y) = -Vjtf{x-y), 



(100) 



V^c ( x — y) being given in the formula (43). Introduce the hermitian bi-local 
bosonic fields $^ (x, y) and the functional integral 

J [£><&] expj-i Jdxj dy<S>£(x,y)U%S (x-y)$<B(y,x)y (101) 

By shifting the functional integration variables, we obtain 

cxp 1 1 J dx J dylp A (x) iP B (y) U%g (x - y) if (y) ^ D (x) } = (102) 

= ~k\ ^ exp {-^ / dx J (x,y)U%S(x-y)$% (y,x)} 



exp 



-rA 



dyip (x) tp B (y) A A (x, y) \ , 



where 



A B A {x,y) = U A 3 g{x-y)$ c D {y,x), (103) 
and rewrite the partition function in the form (78) with the effective action 

S eS [$} = -^ Jdxj dy*i(x,y)U2g(x-y)$Z(y,x) + W[A], (104) 



exp{W[A}} = l + J2 r{n) [ A \ 

71=1 

<„)*,<,)A5(,, 



y)i\ )■ (105) 

Calculations give 

[AI = W (1) \A] = { dx f dy A B A (x, y) S4 (y-x), (106) 



r< 2 > [A] = 1 (V« [A]) 2 + Vt^ 2 > [A] , (107) 

W^[A] = -\j dxj. Jd Vl Jdx 2 J dy 2 A*l(x 1 ,y 1 )Sp i (y 1 -x 2 ) 

A B A l{x 2 ,y 2 )S^\{y 2 -x 1 ), (108) 

r< 3 ) [A] = 1 (V (1) [A]) 3 + TU (1) [A] TU (2) [A] + V^ (3) [A] , (109) 

[A] = lj d Xl J d yi ... j dx 3 J dy 3 A B A \ (x uyi ) S^l ( yi - x 2 ) 

A B A l < 3 - ss) (s 3 , 2/3 ) S£ (» 3 - *i) , (HO) 



The field equation is derived from the variational principle 

SS cS [$] 



<5$g (y,x) 

and has the form 



(111) 



Using the expressions for W^ n ^ [A] , we obtain the explicit equation 
AS(y,x) = Ug A 3 {x-y){S*{y-x) 

- j "dx 2 j dy 2 S^ (y - x 2 ) A% (x 2 ,y 2 ) S£ 3 (y 2 - x) 

+ j ' dx 2 J dy 2 J dx 3 J dy 3 S£ (y-x 2 )A*l (x 2 ,y 2 ) (113) 
S*l (y 2 - x 3 ) A^l (^3, 2/3) S£ 3 (y 3 - x) 

...}. 

Introduce the two-point Green function of quark field in the presence of the quark- 
antiquark pairing 

(y, x) = Si{y-x)- J dx 2 j dy 2 (y - x 2 ) A% (x 2 ,y 2 ) (y 2 - x) 

+ j dx 2 j dy 2 J dx 3 J dy 3 (y - x 2 ) A% (x 2 , y 2 ) S% (y 2 - x 3 ) 

All(x 3 ,y 3 )Sl(y 3 -x) + ... (114) 
It satisfies the Schwinger- Dyson equation 

(y, x) = S£(y-x)- J dz j dw S C B (y - z) Ag (z, w) (w, x) . (115) 

Then the field equation (113) becomes 

A%(y,x) = U8%(x-y)Sj l (y,x). (116) 

Consider the solution of this field equation in the special class of functions de- 
pending only on the difference of the coordinates 

Af {x,y)=A B A {x-y). (117) 

In this case the function (x, y) depends also only on the difference of the coor- 
dinates 



Performing the Fourier transformation 

A B A ( x -y,T-<7) = ±Xy e - '-i-g / rfp e ^ x ^Af (p, £m ), (119) 

P m ( 27r ) J 

e m = (2m + 1) ^, 

and introducing the matrices A (p, s m ) with the elements A^ (p, e m ), we rewrite the 
integral relation (114) and the integral equation (115) in the form of the algebraic 
ones: 

S(p,£ m ) = S(p,e m )-S(p,e m )A(p,e m )S(p,e m ) 

+S (p, e m ) A (p, e m ) S (p, e m ) A (p, e m ) S (p, e m ) ... (120) 

and 

S (p, e m ) = S (p, e m ) - S (p, e m ) A (p, e m ) S (p, e m ) (121) 

or 

1 ' + A(p,e m ). (122) 



S(p,£ TO ) 5(p,e m ) 
The field equation (116) becomes 



(p, e m ) = \ ]T 7^3 / to&CA (P - q> £m ~ e n ) (q, e n ) (123) 

P n ( 27r ) ^ 

where U^a (p — 1' £ ™ -e n ) is the Fourier transform of Uq% (x — y), 

(x-y,r-<7) = jE e ™ mM i/ rf P eiP(y " X) ^(P^»). 



rr-DB , 
U CA 



to m = 2m-. (124) 

Using the field equation (112), we obtain the value of the effective action (104) 

SW [A] 



Seff[*] = W[A]-\JdxJ dyA B A (x, 



y) 



-^jjdxj dyA B A (x,y)S£(y-x) 

+ G ~ I) I dxi I dyi - 1 dX3 J dmA ^ ( Xi 'ri s £ ^- x ^ 

A^ (x 2 ,y 2 ) S£ (y 2 - x 3 ) Af « (x 3 , y 3 ) S£ (y 3 - x x ) 

~(l~2~) J dXl J dyi - J dxi J dy^M { x ^Vi) S bI (Vi ~ x z) 

A% (x 2 ,y 2 )...A*l (x 4 ,y 4 )S£l(y 4 - Xl ) 
+... (125) 



j dx J dyA (x,y)S(y-x) 
J dx 1 J dy 1 ... J dx 3 J dy 3 A (x 1 , y x ) S {y x - x 2 ) 



1 1 

3 ~ 2 

A (x 2 ,y 2 ) S (y 2 - x 3 ) A {x 3 ,y 3 ) S (y 3 - x x ) 

J dxi J dyi... J dx 4 J dy 4 A (x lt yi) S (t/i - x 2 ) 

A (x 2 ,y 2 ) ... A (x 4 , y 4 ) S (y 4 - xi) + ... j , 

where A(x,y) is the matrix with the elements A A (x,y). For the fields in the 
special class (117) the effective action is expressed in terms of the free energy density 
F[x;A] 

S eff m = -3 /dxF[x;Al. (126) 



It follows that 
F[x;A] = F[A] 
- (l- 



+ 



2 

1 I 

3 ~ 2 



| dyA B A 

y*dy y*da;2 J dy 2 J dx 3 J dy 3 A^ (x - j/i) S^ 2 (yi - x 2 ) 
A% - 2/2) S£ (y 2 - x 3 ) Af » (a* - 2/ 3 ) S£ (y 3 - x) 

I " k) I dV I dX2 I J dXi J dy * A * {X _ Vl) S * {Vl ~ X2) 
All - 2/2) - S± (2/3 - x 4 ) All ^ - ^) S£ (J/4 - x) 



+... 
Tri I I 



dyA (x — y) S (y — x) 

+ (1 ~~ 2~) J dy J dx2 J dy2 J dx3 J d V3 A ( x - Vi) s (yi -xi) 

A (x 2 - 2/2) S 1 (2/2 - 2:3) A (2:3 - 2/3) S (2/3 - ar) 
~ (l ~~ 2~) J dy J dx2 J dy 2 --- J dx 4 J dy i A(x-y 1 )S{y 1 -x 2 ) 

A (x 2 - 2/2) ... 5" (2/3 - a*) A (x 4 - 2/4) 5* (2/4 - x) + ...} (127) 
^E^/ dpTr{(l-^jA(p,e m )S(p,e m ) 

~ ~ -1 3 

A(p,e m )S'(p,e m ) 

}• 



'1 1 
+ I 3-2 



+ ... 



I ~~ V A (P' £ ™) S '(P' e ™) 
Summing up the infinite series, we obtain 



> £ m)d£ - -S(p,e m ) 



(128) 



where S ? (p,e m ) is determined by the equation of the form (122) with the replace- 
ment of A (p, e m ) by £A (p, e m ) 



S«(p, e m ) S(p,e m ) 



+ £A ( p , Em ) . 



(129) 



Introduce 4-vector p M with j>4 = e m and denote A^ (p,£ m )by A^ (p) . These order 
parameters have following most general form 

M(P) = S^i^ + ^A^U^ + ilXKtS^ 

+ (7„7 B )f Ajgg („) + {a ^f a aJW^. („) (130) 

The existence of the non-vanishing order parameters with definite transformation 
properties would mean the spontaneous breaking of the corresponding symmetries. 

4 Formation of triquarks 

It is straightforward to generalize the method presented in preceding Sections for 
applying to the problem of the formation of triquarks - the bound states of three 

nnarks We note that the fundamental interaction mechanisms in OC!D Hnstantnn 



induced, gluon exchange, etc..) always lead to some effective (non-local, in general) 
six-fermion interaction between quark fields with the effective interaction action of 
the general form 



Sint = i J dx J dy J dz J du J dv j dwip F (w) tp E (v) tp D (u) 

Vdef {u, v, w; z, y, x) ip A (x) i\) B (y) tp c (z) , 

v def {u,v,w;z,y,x) = -VEef {u,v,w; z,y,x) = -VEEf (u,v,w; z,y,x) = ... 

(131) 

The form-factors VEEf {u,v,w;z, y, x) depend only on the differences of the space- 
time coordinates. The partition function of the many-quark system with this effec- 
tive six-quark interaction equals 



Z = j [Dtp] [Dlf\ exp j- J dxip A (x) L B A i) B (x) 



exp | ^ J dx J dy J dz j du J dv J dwip F (w) ip E (v) ip D (u) 

VEEf («, v, w; z, y, x) ip A {x) 4> B (v) 4>c (*)} ■ (132) 

In order to describe the triquarks we introduce some tri-spinor tri-local field 
^ abc (x, y, z) as well as its conjugate $ (z, y, x) and set 

= j [£>$] [DW\ exp |^ j dx j dy j dz j du j dv j dw 

<S> FED (w, v, u) VEEE (u, v, w; z, y, x) $abc {x, y,z)j. (133) 
By shifting the functional integration variables 



<S>ABc(x,y,z) -» ^ABc(x,y,z) + -^=^ A (x)tP B (y)tP c (z), 

—FED , -FED 1 —F —E , —D 

<P (w,v,u) — > 9 (w, v,u) -\ — -=ip [w)ip (V)1p (U) 

V2 

wc establish the Hubbard-Stratonovich transformation 



exp | i J dx J dy J dz J du J dv J dwip F (10) ip E (v) ip D (u) 
VdEf {u, v, w; z, y, x) ip A (x) tp B (y) i\) c (z)} 
J [£>$] [£>!] exp|-iy dx J dy J dz J du J dv J dw 

<5> FED (w, v, u) VEEf («, v, w; z, y, x) <5> ABC {x, y, z)} 
+lf (z) (y) i> A (x) A ABC {x, y, z)] } , (134) 



where 



A CBA (z,y,x) = J du J dv J dunp F (w)ip E (v)ip D (u)VEEf ( u ,v,w; z,y, x) ,(135) 
A ABC (x,y,z) = [ du [ dv f dwVREE (u,v,w\z,y,x)ib n {u)\b K {v)ib F (w) , 



and transform, after lengthy calculations, the partition function (132) into the form 

QBA 

(23) of a functional integral over the tri-local fields &abc (x, y, z ) and 3> (z, y, x) 
with the effective action 

S eff [$,$] = ~ j dx j dy j dz$ FDE (w,v,u)Vgg£(u,v,w;z,y,x)<f> A Bc(x,y,z) 
+W[A,A\ (136) 
and the functional W [A, A] of the form 

oo 

W [A, A] =J2 w(2n) [ A >^] . ( 137 ) 

n=l 

where W 12 ^ [A, A] is a functional of the n— th order with respect to each type of 

QBA 

tri-local fields Aabc (x, y, z)and A (z, y, x), for example 

W (2) [A, Al = — J dx J dy J dz J du J dv J dw (138) 
A FED (w, v, u) Sp (w - z) S| (v - y) (u - x) A A bc {x, y, z) , 



VF (4) [A, A] = -\j dxij d yi J d Zl ...J du 2 j dv 2 J dw 2 A FlElDl (w^vum) 

^F 2 E 2 D 2 (w2 ^ _ g Bi _ m} gM (ui _ xi) gC 2 ( m _ Z2) 

SeI («2 - 3/2) (W2 - a; 2 ) Aa 2 b 2 c 2 (x 2 ,y2, z 2 ) A Ai b 1 c 1 (xi, 2/1,21) (139) 



From the variational principle we derive the field equation 

]_ ['['[' (5V|^^^^ j~ j 

-A ASC (z,y,;z) = du dv dwV££c (u,v,w;z,y,x)Y^ —fed 

J J J n=1 <5A (w,v,u) 

(140) 

Using explicit expression of W^ 2 ™' [A, A] , we have shown that up to the 6th order 
of the pertubation theory there exists following system of integral equations 

A A Bc(x,y,z) = J du J dv J dw J dx' J dy' J dz'V^Bc (u,v,w;z,y,x) 
G% (w, z') G% (v, y') Gi («, x') A a >b>c> (x' , y' , z') , (141) 

G D A (x, u ) = S A J (x-u) + j dy j dzS B (x - y) Sg (y, z) G° (z, u) , (142) 



(z, w) = J dx J dy J du J dvA FED (w, v, u) Gp (u, x) G E (v, y) A A bc (x, y, z) . 

(143) 

It is easy to verify that (x, u) is the two-point Green function of the quark field 
in the presence of its interaction with the "external" tri-local fields Aabc {x^tz) 

and A (w,v,u). It is determined by the Schwinger-Dyson equation represented 
by the Feynman diagram in Fig. 1. with the self-energy part (143) reprsented by 
the Feynman diagram in Fig. 2. 

The solutions of the system of integral equations (141)-(143) in the class of the 



the anticommuting order parameters of the ground state of the many-quark system 
with the binding of the quarks into the triquarks. This means that in the QCD 
dense quark matter there might exist a phase transition with the anticommuting 
order parameters. Note that there is no condensation of the triquarks, because these 
composite particles are fermions. 
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